Anomalous reaction-transport processes: the dynamics beyond the Mass Action Law 
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In this paper we reconsider the Mass Action Law (MAL) for the anomalous reversible reaction 
A B with diffusion. We provide a mesoscopic description of this reaction when the transitions 
between two states A and B are governed by anomalous (heavy-tailed) waiting-time distributions. 
We derive the set of mesoscopic integro-differential equations for the mean densities of reacting and 
diffusing particles in both states. We show that the effective reaction rate memory kernels in these 
equations and the uniform asymptotic states depend on transport characteristics such as jumping 
rates. This is in contradiction with the classical picture of MAL. We find that transport can even 
induce an extinction of the particles such that the density of particles A or B tends asymptotically 
to zero. We verify analytical results by Monte Carlo simulations and show that the mesoscopic 
densities exhibit a transient growth before decay. 
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The Mass Action Law (MAL) plays a very important 
role in a large number of chemical, biological and physical 
systems [jj, |2|. It states that the rate of an elementary 
reaction is proportional to concentrations of reactants. 
MAL also gives the expression for the equilibrium con- 
stant which is a main characteristic of chemical equilib- 
rium. To illustrate this, consider the reversible reaction 
A ^ B, with (3\ and j3 2 denoting the forward and back- 
ward rate constants. According to MAL, the balance 
equations for the mean concentrations n\ and n 2 of dif- 
fusing particles A and B can be written as 
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where D\ and D 2 are the diffusion coefficients of the par- 
ticles A and B respectively. It follows from (Q]) that the 
uniform equilibrium state (nj°,n§°) obeys the equation 
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where K eq represents the equilibrium constant of the re- 
action process. This constant depends on the thermody- 
namics properties of the system, but is independent from 
the transport parameters D\ and D 2 . The purpose of this 
paper is to reconsider these two fundamental equations 
(fTj) and ([2]) for anomalous reaction and transport. 

Continuous time random walk (CTRW) models have 
been widely used in recent years to gain insights into 
the anomalous transport 0]. The extension of CTRW 
models to reaction-transport phenomena presents mod- 
eling challenges, because of the difficulty of taking into 
account chemical reactions within anomalous transport. 
Recently several authors have explored the reaction- 
transport models in which the standard diffusion is re- 
placed by an anomalous (subdiffusive) transport 0, [U, d, 
0, H, [tj O) EH- I* nas been shown that the evolution 
equations for density of particles are drastically differ- 
ent from the standard reaction-diffusion equations. For 



example, the transport and reaction terms are not sep- 
arable as it happens in the classical case (JI]). Instead, 
one finds that the transport term becomes dependent on 
the reaction constants [3\ or (3 2 [1, Q . The master equa- 
tion for the mean densityof one of reactants may include 
crossed transport term p], |Tl[ . This is a consequence of 
the non-Markovian nature of subdiffusion. 

In previous works, however, reaction was always in- 
troduced phenomenologically following the principles of 
classical reaction kinetics. The idea of this paper is to 
consider both the reversible reaction A B and sub- 
diffusive transport from a probabilistic point view. It is 
well known that the classical kinetics in |T]) corresponds 
to Markovian transition of particles from one state to an- 
other. Our aim is to take into account anomalous (non- 
Markovian) transitions of particles from the state A to 
B and backward and find out how the transport process 
and reactions are coupled. In what follows we will show 
that fundamental constant K eq becomes dependent on 
transport characteristics which is in contradiction to the 
classical picture of MAL. This is due to anomalous nature 
of transitions A <^ B for which the waiting time distribu- 
tion exhibits the power law decay with the infinite mean 
waiting time. Let us mention that the situations that are 
outside the scope of the MAL have already been reported 
for diffusion-limited reactions with long-range interac- 
tions in space [HI • The reversible reaction At± B can be 
interpreted as a switching between two states A and B. 
This topic has attracted a great interest recently because 
the switching process can be non-Markovian. Examples 
include two-state ion channel gating [l~5l ]. stochastic res- 
onance [IH , quantum dots [iji [T|| , etc. For anomalous 
switching process without transport the mean residence 
time of the particles in each state is divergent [l7], Ell • 
As a result, the density of particles in one state tends to 
zero in the limit t — > oo which means the extinction of 
one of the states. Here we show that the transport pro- 
cess can drastically modify extinction/survival dynamics 
for anomalous transitions. One of the motivations for our 
study is the experimental data for a malignant brain can- 
cer that exhibits migration-proliferation dichotomy [l|. 



The motility (transport) of cells and phenotype transi- 
tions A <=^ B (prolifcration^migration) can be anoma- 
lous simultaneously Another possible application of 
our model is the isomerization reaction for which macro- 
molecule in two interconvertible states migrate with dif- 
ferent electrophoretic mobilities fljjj ]. 

In this Letter we consider the following stochastic 
model for the transport and reversible reaction A ^ B. 
The particles of type A and B randomly move along one- 
dimensional space and switch between the states A and 
B. This random walk with switching can be described 
by four sequences of mutually independent random vari- 
ables. Two sequences {YiY 2: ...} and {Zi^Z^,...} describe 
the waiting times between jumps for particles in the state 
A and B correspondingly. We assume that these ran- 
dom variables are identically distributed with probabil- 
ity density function (pdf) (fi{t) for particles A and pdf 
(f 2 (t) for particles B. Two other sequences {Ux t U 2i ...} 
and {Wj. W2 .} describe the waiting times for random 
transitions: A — * B and B — > A respectively. The ran- 
dom variables UiU 2t ... and Wi ) W2,...are identically dis- 
tributed with pdf's ipi 2 (t) and ip2i(t)- If we place the 
particle at the position x at time in state B 7 and if the 
random time Z\ for a jump is less than random time W\ 
for reaction B — > A, then the random jump happens at 
time Z\ . However if W\ < Z\ , then the transition B — > A 
occurs at time W\. In other words an event (reaction or 
jump) happens at time min (W\, Z\) . For example, there 
is a jump at the position X\ at time Z\ < W\, then a 
second jump of length X2 after a further time Z2 < W2, 
a switching to state A after time W3 < Z%, a transition 
back to B after time U4 < Y4, and so on (see Fig.l). 
So jumps and transitions A ^ B are not independent 
random events as in [7], [ll[ . 

We express the density of particles Uj(x, t) at position 
x at time t in terms of the initial distribution of particles 
n1(x) and density of particles ji(x, t — t) that arrive at 
the same position x at previous time t — r. The balance 
equations for the mean densities n\ and ri2 of particles A 
and B are 

ni(x,t) = nS(s)$i(t)*ia(t) + 



ji(x,t - t)$i(t)#i 2 (t)<2t, 



n 2 (x,t) = nSj(aj)$ 2 (t)*2i(t) + 



ia(a:,t-r)$ 2 (r)*2i(r)dr, (3) 



where = (/?i(r)dr and *y(t) = ipij(r)dT 

are the corresponding survival probabilities for ^ (i) and 
ipij(t). For example, $i(f) is the probability that a parti- 
cle in the state A does not jump until time t, and ^i 2 (i) 
is the probability that a particle in the state A does not 
switch to B until time t. The first term in the RHS of the 
Eq. (|3|) represents the contribution from the initial den- 
sity of particles that have neither jumped nor switched 
until time t. The density ji(x,t) describes how the par- 
ticles arrive at point x at time t as a result of the trans- 
port and switching processes. Equations for the density 
ji(x, t) of particles A and the density j 2 (x, t) of particles 
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Fig. 1: A single realization of the random walk with switching 
A ^ B. Initially the particle is assumed to be in the state B 
and at position x = 0. 



B can be written as 
j x {x,t) = n° 2 (x)$ 2 (t)il>2i(t)+ 
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The hrst equation is the conservation law for particles 
of type A at time t at position x. The first term in the 
RHS of the equation accounts for the initial distribution 
712(0;) of particles in state B that switch to A at time 
t, provided they do not jump up to time t. The second 
term is the contribution from the initial density n\(x) of 
particles in state A that jump to x from x — z at time 
t having not switched until t. The random jump length 
z is determined by the dispersal kernel pi(z). The third 
term represents the contribution from those particles that 
switch from state B to state A after a waiting time t, un- 
der the condition that they do not jump during that time. 
Finally, the fourth term corresponds to the contribution 
of particles in state A that arrive at x — z at t — r do not 
switch to B during time r. 

The set of linear equations ([3]),([4]) can be solved by 
using the Laplace-Fourier transforms: (x,t) — ► (q,s). 
One can obtain two equations: sni(q,s) — n\(q) = 
h (s) (pi (q) - l)ni(<?, s) - ai (s) n\{q, s) + a 2 (s) n 2 (q, s) 
and sn 2 (q,s) - n%(q) = k 2 (s)(p 2 (q) - l)n 2 (q,s) + 
ai (s) ni(q, s) — 0,2 (s) ni(q, s). Here we introduce 



ki (s) 
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with the notation [/] = f (s). Taking the inverse 
Laplace and Fourier transforms we obtain the following 



master equations 
dm '■' ' x 



-gj- = I I ni(t-r,x- z)ki(T)(pi(z) -l)dzdr+ 
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Oi(r)ni(t — r, x)dr + / a 2 (T)n 2 (t — T,x)dr 



dn f* f 00 

= y J n 2 {t - t,x - z)k 2 (T)(p 2 (z) - l)dzdr+ 



+ / ai(r)ni(t — r,x)dr — / a 2 (r)n 2 (t - r, x)dr, 
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where fcj(i) and aiit) are the inverse Laplace transforms 
of ki(s), <n(s) defined in (0. The most interesting 
feature of the system © is that effective reaction rate 
memory kernels a\ (t) and a 2 (t) depend on the transport 
through the survival probabilities &i(t) and & 2 (t), while 
the transport memory kernels ki (t) and k 2 (t) depend on 
statistical characteristics of reactions such as tpij . If the 
random waiting times for switching and jumping are ex- 
ponentially distributed: ipi — \e~ Xit , ipij — f3ie~^ it , 
then these dependencies cease to exist. As a result the 
transport and reaction terms are separable as in the clas- 
sical case ([1]). For example, if we use the diffusive ap- 
proximation for transport pi (q) ~ 1 — af q 2 , then the 
system © can be written as classical reaction-diffusion 
equations (JTJ) , with Di = Xiof. Similarly, for a Marko- 
vian switching process with subdiffusive transport, we 
could recover from ^ the model for cancer spreading 
studied in Q . If the waiting time pdf ip\ 2 has a gamma 
distribution as tpi 2 = j3\te~^ xt and ipi — Aie~ Al *, then 
ai (s) = 01 (2/3i + Ai + s) . So the reaction rate mem- 
ory kernels are 



i = 1,2. 



(7) 



This formula shows that the effective reaction rate kernels 
depend on the rate of jumps A^. Now, let us find the 
uniform stationary states corresponding to ^ under the 



condition 



1. From the limit q — > one finds 



Pi (q) — 1, so we obtain the asymptotic state as s 
a 2 (s) a x (s) 
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1 1 " \ai (s) + a 2 (s) ' ai (s) + a 2 (s) y 

The main feature of this asymptotic state is that in gen- 
eral it depends on the characteristics of the transport 
process which is in contradiction with the Mass Action 
Law @. This follows from the fact that the survival 
function <i>i appears in the definition of cij. This happens 
for any situation except when the switching process is 
Markovian for which lim s _>o en — Pi . 

Assume now that the reaction process is governed by 
a power-law decaying distribution of waiting times. We 
use the standard approximation [ipi 1 1 ' 1 



i - fa 1 *) 



with jij < 1 as s — *■ 0. On the contrary, for the transport 
we consider the Markovian case: <fi(t) — A;e~ Ait . Then, 
we find that the state (JSJ) turns into 
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Fig. 2: Time evolution of the density ni obtained from Monte 
Carlo simulations for different values of the parameters Ai and 
A 2 ; o\ = a\ = 1, /3i = p 2 = 1, 7 i2 = 0.25, 721 = 0.5. 



where the explicit dependencies of asymptotic states on 
the transport parameters Ai and A2 are evident. In fact, 
the ratio of the two uniform densities in the limit t — > 00 
is 



^X— 721 yJ72i 
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= K, 



(10) 



where K becomes dependent on the transport parame- 
ters Ai and X 2 (see © for comparison). Note, however, 
this new constant K cannot be considered an equilib- 
rium constant since in the non-Markovian case consid- 
ered here a thermodynamic equilibrium state cannot be 
defined. Let us assume (3\ < (3 2 ; then the MAL ^ pre- 
dicts nj° > nf. On the contrary, it is clear from (fTU|) 
that one could choose the rates of jumping Ai and X 2 so 
that the inequality can be inverted to nf < n™. We 
obtain even more dramatic results if we take the limit 
Ai — > or A 2 — > 0. If, for example, we consider the limit 
Ai — > (A 2 7^ 0), then one can observe the extinction 
of particles in state B and survival of particles in state 

A, that is, m — * 1 and n 2 — > as t — > 00 (see ©). 
So we find from our model that transport process can 
induce a survival/extinction of one of the two densities 
for anomalous reactions. To validate this phenomenon 
we have performed the direct Monte Carlo simulations 
of two-states random walks. The results are illustrated 
in Fig. 2 where one can see that if Ai and A2 = 0, 
then we might observe the temporal growth of n± before 
the final decay to zero (solid line). However, if we put 
Ai = A2 = (dashed line), then the limit for the den- 
sity of particles in state A is completely different, that is, 
m — > 1 as t — > 00. 

The result m — > 1 as Ai — > imply that if the parti- 
cles do not move in one of the states, they survive. This 
idea of 'staying quiet helps you to survive ' can be under- 
stood from the interplay between the waiting times for 
reactions and jumps. According to our derivation, the 
reaction process, say the transition from state A to state 

B, is actually governed by the density $1^12 (the parti- 
cles react only if they have not jumped before, as can be 
seen from (gj). We can refer to this function as the effec- 
tive waiting time pdf. For an anomalous switching pro- 
cess with Markovian transport the asymptotic behavior 

(9) of the effective pdf reads $1^12 ~ i _1 " 7l2 e" Al *. Then, 



2 



+ 1 



Extinction of B 
(n- = 1; n 2 - = 0) 



Coexistence 
of A and B 



Extinction of A 
(a' = 0; n 2 ~ = 1) 



1 



Fig. 3: Phase diagram of the extinction/survival regions for 
anomalous switching and transport in terms of the exponents 
ji and jij. 



the mean waiting time is finite, and for this reason the 
system reaches a stationary state, given by ©. How- 
ever, in the limit Ai — > the effective mean waiting time 
diverges, which makes the reaction A — ► B much slower 
than the backward reaction B — > A, so the particles tend 
to get trapped in the state A. For this reason we obtain 
n\ — > 1 and ri2 — > in the long-time limit. 

For anomalous transport and anomalous reaction 
(switching), we have ipij(s) ~ 1 — (/3~ 1 s) 7lJ and (fi(s) ~ 
1 — (A^ 1 s) 7 *, with 7< < 1 and 7^ < 1. ft is helpful again 
to use the idea of an effective waiting time distribution 
for reaction: Qiipij ~ w - ^. So that, the effective 
mean waiting time can be finite, provided that the con- 
dition 7i -(- 7y > 1 is satisfied. Fig. 3 shows the "phase 
diagram" for asymptotic states nj° and depending of 
the values of 71 -I- 712 and 72 + 721- If both of them are 
larger than one, then the transitions 1 — > 2 and 2 — > f 
are governed by finite mean waiting times, so a coex- 
istence of two states is possible. In other regions, the 
divergences of the mean waiting times make the parti- 
cles get trapped in the state where the switching process 



is slower. So that, one of the states become extinct in 
the asymptotic regime. These results can be explained 
by a coupled renewal property assumed in our model. If 
the "internal" waiting time of the particles starts from 
zero after each event (reaction or jump), then we have a 
competition between both processes to be the first to oc- 
cur, and so coupled effects emerge. This coupled renewal 
property is opposite to additive renewal property when 
the random walk in space is completely independent of 
the reaction process (see, for example, [7|). 

To sum up, we have presented a non-Markov model 
for the reversible reaction A <=^ B and studied the in- 
terplay between anomalous transport and anomalous re- 
action process implemented in a probabilistic way. So 
we have been able to explore those situations that are 
beyond Mass Action Law. We have derived mesoscopic 
integro-differential equations for the mean densities of 
particles in states A and B when the transitions between 
two states A and B and jumps in space are governed 
by heavy-tailed waiting-time distributions. It has been 
shown that the anomalous properties of the reversible 
reaction yield the appearance of unusual properties such 
as dependence of asymptotic states on transport and the 
transient growth of densities before decay. We have found 
that the transport process can modify completely the uni- 
form stationary regimes. In particular, it can induce the 
survival/extinction of one of the states. These results 
have been validated by direct Monte Carlo simulations 
of two-state random walk. 
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